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ABSTRACT 

We investigate how the scaling relations between central black hole mass and host galaxy 
properties (velocity dispersion, bulge stellar mass and bulge luminosity) depend on the large 
scale environment. For each of a sample of 69 galaxies with dynamical black hole measure- 
ments we compile four environmental measures (nearest neighbor distance, fixed aperture 
number density, total halo mass, and central/satellite). We find that central and satellite galax- 
ies follow distinctly separate scalings in each of the three relations we have examined. The 
JVI, - cr relation of central galaxies is significantly steeper { j3 = 6.39 ± 0.50) than that of satel- 
lite galaxies (/? = 4.78 ± 0.51), but has a similar intercept. This behavior remains even after 
restricting to a sample of only early type galaxies or after removing the 8 brightest cluster 
galaxies. The M, - a relation shows more modest differences when splitting the sample based 
on the other environmental indicators, suggesting that they are driven by the underlying satel- 
lite/central fractions. Separate relations for centrals and satellites are also seen in the power 
law scaling between black hole mass and bulge stellar mass or bulge luminosity. We suggest 
that gas rich, low mass galaxies undergo a period of rapid black hole growth in the process of 
becoming satellites. If central galaxies on the current M, - a relation are representative pro- 
genitors of the satellite population, the observations imply that a cr = 120 kms^^ galaxy must 
nearly triple its central black hole mass. The elevated black hole masses of massive central 
galaxies are then a natural consequence of the accretion of satellites. 
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> 1 INTRODUCTION 



The correlation between the mass of a central blacli hole (BH) 
and the properties of its host galaxy suggests black holes play 
a major role in galaxy evolution and has stimulated a wide 
range of theoretical and observational studies. It is now clear 
that the mass of the black hole forms a tight relation with 



the velocity dispersion of its galaxy (jpe rrarese & MerritI 
Gebhardt et al. 2000; Tr emaine et"al1 [200 2; Graham et al 



200C 



2011 



McConnell&Ma 2013), the mass of its bulp dMagorrian et al.l 
1998; Marconi & Hunt 20031; iHiiring & Rixll2004l). the luminosity 
of its bulge I Kormendv & Richstondl 199^ Kormendv & Gebhardtl 



l200ll ;|Graham"2007'; 'Graham & Scott"2013') and the Sersic profile 
of the bulge (Graham et al. 2001; Graham & Driver 2007). It has 
been postulated that these tight relations point to a self-regulation 
between the growth of the galaxy and its black hole. Feedback from 
an active galactic nucleus (AGN) is a plausible mechanism for this 
self-regulation and may explain the quenching of massive galax- 
ies and the properties of groups and clusters (Silk&Rees 1998; 
Kauffmann & Haehnelt 200G; King 2003; Springel et al. 2005; 
Croton et al. | i2006i; iBower et al.i i2006.; .Hopkins et al., ,2008b 



Croton et al.i i2006i; iBower et al.i ,2006.; .Hopkins et al., ,2008biiai ; 
Booth & Schavell2009| - |McCarthv et al.l201(]| ; (Dubois et alj2012ah . 
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Nonetheless, many questions remain about how th e AGN man- 
ages to accrete gas ( Cattaneo & Teyssier 2007; Hopkins et akl 
l2006bl ; |Pizzolato & Sokeiii2010l ; lDubois et al. 2012b) and how the 
AGN outflow couples t o the surroundi ng gas (.Murray et al.l 



Ciotti & Ostrike^ (20071; los triker et all i2010l: 



Novak et 



Faucher-Giguere & Ouataerlil2012i;iSilk & Nusseij|2010 l) 



2005 



2012 



A powerful way to test models of black hole growth is by 
examining how their scaling relations change in different environ- 
ments. It has been known for some time that satellite galaxies, or 
more generally galaxies in dense e nvironments, tend to be pas 



sive, spheroidal and lacking in gas (lDressleij|l9 80; Havn es et al 
Il984l; Postman & Gelle3ll984l; iBalogh et al.lll999; .Blanton et al 
2005; Poggianti et all |2006|; Iwilman et al.l |2009|; iBamford et al 
120091 ; |Pengetal.ll2O10l ; iMcGee et alj 1201 ll ; IWetzel et alj |2012|) . 



It is now clear that galaxies being accreted into dense environ- 
ments must halt their star formation and alter their morphol- 
ogy. Although the exact physical mechanism responsible is not 
known, it seems to be related to the gas stripping that a satel- 
lite undergoes as it encounters the hot medium within groups and 
clusters (Gunn & Go tt 197 2; Abadi et a l. 1999; Ouilis et al. 200(i|; 
[McCarthy etaLll2008l ; lBahe et al II2OI2I) . The morphological trans- 
formation associated with becoming a satellite galaxy is especially 
important for the black hole scaling relations because it may pro- 
vide a new channel to the formation of bulges and early type galax- 
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ies. Unfortunately, how satellite galaxies transform from spirals 
to early types is a matter of some disagreement in the literature. 
It has been suggested that the bulges of satellite early types are 
too large to come from s piral galaxies w i thout some additional 
bulge growth {c hristlein & Zabludoffll2004lBekld & Couchll201 ll : 
iLackner & Gunni,2013i) . However, large samples of resolved color 
profiles are compatible with a gentle fading of disks and the reveal- 
ing o f their prexisting bulges dWeinmarm et al.ll2009l : [Hudson et alj 
I2OIQ) . Either way, this unique route of galaxy evolution occurs only 
in satellites and thus is a powerful testbed for the nature of black 
hole scaling relations. 

The distinction between a central galaxy and satellite galaxy 
may be most important for understanding galaxy evolution but it 
is often difficult to measure definitively. Further, there are reasons 
to believe as many as three types of 'environment' other than cen- 
tral/satellite might be important for black hole scaling relations. 
Firstly, the processes of harassment and merging rely largely on in- 
teractions with other galaxies, which implies that the likelihood for 
interaction, or the distance to an Nth nearest galaxy, is important. 
Secondly, there is evidence that dark matter halos form sooner in 
large scale overdensities on the order of ~ 5-10 Mpc jLi et alj2008l : 
[Oao et al. 2005). If this growth is directly linked to the formation 
of galaxies, then it is possible that galaxies in large scale overdensi- 
ties have a 'head-start' in their development. Finally, the total mass 
of a galaxy's host halo, regardless of the galaxy' s location, may in- 
fluence how the ga s is accreted onto the galaxy jSimha et al.l2009l : 
IZavala et alj2012h . 

In this paper, we calculate several environmental measure- 
ments for e ach galaxy with a direct black hole measurement as 
compiled bv lMcConnell & Mj j2013l) . Each of these environmen- 
tal measurements is designed to probe a potentially important envi- 
ronment. We will use these measurements to divide the sample of 
black holes into subgroups and examine the resulting scaling rela- 
tions. In ij2]we discuss the black hole and galaxy property measure- 
ments we use and how we derive the environmental measures. In Sj3] 
we show the strong effect environment has on the black hole scaling 
relations and in !j4]we discuss how this could arise. We conclude in 
35] In this paper, we adopt a ACDM cosmology with the parame- 
ters; f2m = 0.3, Q.K = 0.7 and h = Ho/(100 kms"^ Mpc"^) = 
0.70. 



2 DATA AND METHODS 

We plan to examine the scaling relations between black holes and 
their internal galaxy properties with the aim of determining the 
effect of their external, or large scale, environment on such rela- 
tions. The large scale environment profoundly affects many galaxy 
properties, such as star formation rate and morphology. Indeed, it 
is this effect which makes examining the black hole scaling rela- 
tions in different environments a worthwhile endeavor However, 
for this reason we must restrict our sample of black hole masses to 
those which are determined as directly as possible, that is, those 
derived from stellar or gas dynamics. Several indirect measure- 
ments, such as those based on the width of emission lines, ex- 
ist (McLure & Dunlop 2004; Shen et al. 2008), but these are cali- 
brated on the BH scaling relations and could also be systematically 
affected by environmental conditions. Thus, we will restrict our dis- 
cussion to low redshift scaling relations where the black hole mass 
has been more directly measured. 



2.1 Black holes and their host galaxy properties 

We take as ou r start ing point the recent compilation of 
iMcConnell & M j ilOlj) . These authors have presented critical 
evaluations of nearby black hole and host galaxy properties accom- 
panied by judicial pruning based on the quality of the measure- 
ments. As a base, we will use the black hole mass, velocity disper- 
sion, V band bulge luminosity and bulge stellar mass of their full 
sample of 72 galaxies. While each galaxy has a measurement of the 
black hole mass and velocity dispersion, only 35 of the early type 
galaxies have bulge stellar mass measurements and 44 have bulge 
luminosities. 

Of this base sample of 72, as discussed in H2.2\ we were not 
able to make environmental measures for two galaxies in the sam- 
ple (the Milky Way and Circinus). Another galaxy, M32, is of sig- 
nificantly lower velocity dispersion, luminosity and mass than the 
rest of the sample and could bias the scaling relations when the 
sample is split in two. Therefore, we remove all three galaxies from 
further analysis, leaving a reduced sample of 69 galaxies. 

2.2 Environmental measures 

In an attempt to probe each of the different environments dis- 
cussed earlier, we will use 4 environmental indicators. Each of 
these indicato rs is either calculate d from the 2MASS Redshift Sur- 
vey (2MRS) jHuchra et alj|20l3) or drawn from previously pub- 
lished analysis by the 2MRS team. The 2MRS was designed to 
obtai n a redshift for all galaxies with Two Micron All Sky Sur- 
vey i Skrutskie et al.ll2006t. 2MASS) Ks magnitudes brighter than 
11.75 and greater than 5° from the galactic plane. As mentioned, 
two galaxies from the original sample of 72 galaxies with good a 
and M. measurements are not in the 2MRS: the Milky Way, for 
an obvious reason, and Circinus, which is less than 5° from the 
galactic plane. Otherwise, 2MRS provides a nearly uniform sam- 
pling of the galaxies in the local universe, and thus is ideal for our 
purposes. While not a strict requirement for our puiposes, the Ks 
selection provides for a nearly stellar mass selected sample. This 
is particularly useful when measuring environments, where the star 
formation histories and thus the optical luminosities systematically 
vary as a function of environment. The distribution of redshift and 
Ks magnitude of the full 2MRS and of galaxies with good M, mea- 
surements is shown in Figure[T] 

When calculating any environmental measure from a flux lim- 
ited survey, care has to be taken to adjust for the varying intrinsic 
luminosity limit as a function of redshift. One common method is 
to make a volume limited sample by restricting the tracer galaxies 
to have absolute magnitude brighter than the flux limit at the maxi- 
mum redshift of the sample. The most distant galaxy in our sample 
is the Brightest Cluster Galaxy (BCG) of Abell 1836 at a distance 
157.5 Mpc, corresponding to a absolute magnitude Ks limit of - 
24.3. This limit is shown as a red line in Figure [T] This method 
is used for the first two environmental measures. Alternatively, the 
magnitude limit can be directly modeled by using a sliding esti- 
mate of the typical inter-galaxy distance. This second method was 
implemented during the creation of the group catalog that we use. 

2.2.1 Nth nearest neighbor distance 

The projected distance to the nearest Nth galaxy nei ghbor has been 
used extensively a s an environmenta l indicator (eg., lDressleJl98"ol; 
iBaldrv et al ] |2006l ; |Penge"tai]|2010l) and is not strongly dependent 
on the choice of parameters. To calculate this quantity, we define a 
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Figure 1. Magnitude vs redshift of the 2MASS Redshift Survey (2MRS) 
and the galaxies with dyanamical black hole measurements. The blue line 
shows the magnitude limit used for the tracer galaxies in the calculation of 
the background environmental measures (4th nearest neighbor distance and 
number density). 



subsample of tracer galaxies from the 2MRS which have a magni- 
tude brighter than the Ks limit of -24.3 and a line of sight velocity 
of ± u from the target galaxy. The projected distance to the TVth 
nearest galaxy in this tracer sample is the quantity of interest. In 
this paper, we set v = 500 kms~^ and take A*' = 4. The distribution 
of the resulting distances is shown in the top panel of Figure|2] This 
is a direct probe of the galaxy density, but is difficult to physically 
interpret as the effective scale of the measurement varies. 



2.2.2 Fixed aperture number density 

A complement to the A^th nearest neighbor distance is a measure of 
the number of galaxies within a fixed aperture cylind rical volume. 
This quantity has also been used in many studies (eg. , ICroton et al.l 
I2OO5I : iBlanton & Berfin j I200I Iwilman et al.l I2OI0I) and for our 
purposes is ideal to measure uniformly the very large scale environ- 
ment of a galaxy. While a typical galaxy cluster may have a virial 
radius of ~ 2 Mpc, there are overdensities on much larger scales (~ 
10 Mpc) which our other indicators do not effectively probe. The 
calculation of fixed aperture number density requires the definition 
of a radius, 7?, and line of sight velocity, ± v, to create a cylin- 
der centered on each galaxy. With these definitions, the number of 
tracer galaxies in each cylinder is counted. In an effort to probe the 
large scales, we define ii as 5 Mpc and v as 1000 km/s. The result- 
ing distribution of number densities is shown in the bottom panel 
of Figure |2] 



2.2.3 Group membership 

The 2MRS team has extensively studied the large scale structu re in 
their survey, presenting both gr oup catalogs ( Crook et al .l2007h and 
a reconstructed density field (iErdogdu et al. 20061) . We make use 
of the high density contrast (H DC) group catalog o f ICrook et all 
( I2OO7I) , subsequently updated in ICrook et aU ( |2008|). This catalog 
uses the well known friends-of-friends method ( Huchra & Gelleri 
Il982h to link galaxies into groups and clusters. In Crook et al., 
the projected separation linking length, D, is scaled with redshift 
to account for the flux limited sample and thus the varying inter- 
galaxy distance. In the middle of the redshift range, cz ~ 5000 



km/s, the linking length is 1 Mpc. The line of sight velocity dif- 
ference is not scale d with redshift and is set to be equal to 350 km/s. 
ICrook et"H] ( 1200 7') show that this choice of linking parameters cor- 
responds to a density contrast of 5p/ p = 80. They also present a 
lower density contrast group catalog (LDC) with larger linking pa- 
rameters. However, this more generous catalog tends to combine 
'known' galaxy clusters, such as Abell 426 and Abell 347, into sin- 
gle objects. In this paper we will use the HDC catalog, but we have 
also done the full analysis using the LDC catalog with no signifi- 
cant change in the results. As an example of how the group catalog 
correlates with our previous measurements, the red line in Figure|2] 
represents the distribution of environmental measures for galaxies 
in halos with total virial masses greater than 1O^'''^M0. Clearly, 
galaxies in these dense environments are equally well described by 
any of our measures so far. The 2MRS group catalog is complete 
only to cz < 10000 kms^^, which leaves three galaxies beyond the 
limits of the catalog (NGC6086, NGC6264 and the BCG of Abell 
1836). For the sake of consistency, we remove these galaxies from 
consideration in this section only. 20 galaxies of the remaining sam- 
ple of 66 galaxies are not linked to groups in the lCrook et al.|j2007h 
catalogue. For these galaxies we will assume that their halo masses 
are below the limits of detection, and are thus in low mass halos. 
These galaxies will all be in the lower 50% samples of group mass. 

2.2.4 Central or satellite distinction 

As already stated, the location of a galaxy within a dark matter halo 
plays a key role in its evolution. As a complement to our previous 
measures, which did not take this into account, we will also calcu- 
late whether a gal axy is a central or satellite. For galaxies within 
ICrook et al.l ( |2007|) , we will assume that the most luminous galaxy 
(Ks band) in each group is the central galaxy. While this definition 
makes no link to the location of a galaxy, it is a common way to 
divide between centrals and satellites. A key reason for this is the 
difficulty of measuring the true center of the group of galaxies from 
their member positions alone. While this may be trivial for a galaxy 
cluster with a thousand members, it is essentially meaningless for 
a group of 3 members. Nonetheless, this definition is reasonable, 
because it picks out the known BCGs of Coma, Abell 3565, Abell 
1367, Abefl 2162 and Abefl 2666. However, NGC1399, which is 
commonly thought to be the central galaxy in the Fornax cluster, is 
actually fainter (Ks = 6.4) than a galaxy near the edge of the cluster 
(NGC1316; Ks = 5.9). Similarly, M87 (Ks = 5.9) is dynamically at 
the center of the Virgo cluster, but M49 (Ks = 5.5), a galaxy at the 
center of an infalling subgroup, is actually the brighter galaxy. For 
consistency with groups which do not have well measured centers, 
we do not change these. Thus, it is important to remember that our 
definition of a central is based on luminosity alone. 

For galaxies which are not members of groups, the determina- 
tion of their central/satellite status is more difficult. While the temp- 
tation is to claim that all non-grouped galaxies are centrals, this is 
clearly not the case. We will assume that a galaxy is a central if it 
is the brightest Ks galaxy within a given radius and velocity range. 
One option for the radius to use would be to scale by the virial 
radius of a typical halo hosting a galaxy of the given magnitude. 
However, while the virial radius of a typical low mass galaxy is on 
the order of a few hundred kpcs, a galaxy cluster may be 2 Mpc. 
So, a galaxy may appear to be a central even while bein g embedded 
in a m assive cluster with a much more luminous BCG. lGeha et al.l 
l l2012h have shown that low mass galaxies have non-zero passive 
fractions until they are more than 1.5 Mpc from a nearby massive 
galaxy. This suggests that a reasonable value for the radius would 
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Figure 2. The cumulative distributions of environmental indicators, (top) 
Shown for the 4th nearest neighbor distance as defined in i|2.2.1| Dense 
environments have low nearest neighbor distance, {bottom) Shown for the 
10 Mpc number density as defined in i|2.2.2| In this case, dense environ- 
ments are characterized by high number density. In both cases, the black 
lines are all 69 galaxies in our BH sample, while the red lines show only the 
distribution of the galaxies in galaxy groups with lO^^'^MQand greater. 

be 1.5 Mpc. For the velocity difference we are motivated by the 
350 kms~^ velocity difference used in the creation of the group 
catalog. Therefore, we call all non-grouped galaxies which have a 
more luminous Ks band galaxy less than 1.5 Mpc and ± 350 km/s 
away a satellite galaxy. This adds an additional 6 satellites from the 
population of 23 non-grouped galaxies. 



2.3 Line fitting 

We parameterizing the power law scaling relati ons between black 
hole mass and other parameters, following iMcConnell & Mai 
( l2013h . as: 

logioAf. = a + /31ogioX, (1) 

where X is either cr/200kms~\ Li^/IO^L© or Mb„igs/1O"M0 
and M, is measured in units of M©. Thus, the slope and intercept 
of each scaling relation is given by /? and a respectively. 

We fit this e q uation using the linear regression Bayesian es- 
timator of iKellvl boOTi) . This algorithm accounts for measure- 
ment error on both variables and also prod uces an estimate of 
the intrinsic scatter, e. lMcConnell & Mii (2013j) found this method 

i ;ave consistent result s with the least squares e stimator HPFITEX Y 
Williams et ZII2OI0I) which builds on MPFIT ( lMarkwarddl2009h . 
There has been much discussion in the literature about how to fit 
the power law scaling relations. While this is critical for examining 
the absolute values of the scaling relations, our goal is to examine 
the differential evolution of subdivided samples and is thus slightly 
less dependent on the details of the regression scheme. We perform 
the linear regressions in the forward direction, but we note that the 



qualtitive trends we see are similar if the linear regression is done 
in the inverse direction. 



3 RESULTS 

We will now examine the results of performing linear regressions 
on each of the three primary black hole scaling relations while di- 
viding into several subsamples based on the calculated environmen- 
tal indicators. 



3.1 M,-cr relation 

First, we look at the full sample of galaxies with which we were 
able to make good environmental measurements, a total sample of 
69 galaxies. We find a = 8.33 ± 0.05 and P = 5.62 ± 0.34 with 
a scat ter, e, of 0.39 ± 0.04. As a comparison, iMcCoimell & Mai 
( l2013h find very similar results of (q, /?, e) = (8.32 ± 0.06, 5.58 ± 
0.34, 0.42 ± 0.04) using the same fitting method, but with the full 
sample of 72 galaxies. Our removal of 3 galaxies has not affected 
the bulk properties of the scaling relation, and we confirm that we 
obtain exactly the same results as iMcConnell &Mj j2013l) for the 
full sample. 

We sort the 69 galaxies by each of 4th nearest neighbor, 5 Mpc 
number density and total halo mass, and perform a linear regres- 
sion using the lowest 25%, the lowest 50%, the highest 50% and 
the highest 25% of each indicator. While we call the percentiles by 
these names, in detail, depending on the number of galaxies, the 
exact percentiles vary slightly throughout the sample. The full re- 
sults, including the intercept, slope and scatter, as well as the num- 
ber of galaxies in each subsample, are presented in Table [Tj The 
top three panels of Figure [3] show the results of fitting the scaling 
relations to the highest 50% and lowest 50% of each of 4th nearest 
neighbor, 5 Mpc number density and group mass. In each case, the 
intercept, q, of the relation in the dense environment (ie, shortest 
4th nearest neighbour distance, highest number density, and largest 
halo mass) is systematically higher than in the more isolated envi- 
ronments. The largest difference occurs in the halo mass division, 
where a is 8.21 ± 0.08 in low mass halos and 8.45 ± 0.08 in high 
mass halos. There is also a tendency for the denser environments 
to have shallower slopes, but because of the large uncertainties in 
determining slope this is not statistically signficant. For instance, in 
more massive halos, the slope is /3 = 5.56 ± 0.63, and in smaller 
halos it is 5.28 ± 0.46. The intrinsic scatter in each subsample is in 
the range e = 0.35 - 0.54 without a significant systematic trend with 
environment. 

In the bottom panel of Figure [3] we show the results of sep- 
arately fitting the satellite galaxies and the central galaxies. Inter- 
estingly, here the results are different from the other environmental 
indicators. The central and satellite galaxies have similar intercepts 
a (8.28 ± 0.07 vs 8.36 ± 0.08) but different slopes. In the sample 
of central galaxies /3 = 6.39 ± 0.50 while in satellites /3 = 4.78 ± 
0.51. The departure from the trends of the other indicators makes 
sense. At high velocity dispersion, central galaxies will be in mas- 
sive halos, while at low velocity dispersion, satellites will be in 
more massive halos. The steep slope of the central galaxy relation 
leads to a vertically offset relation when dividing only be group 
mass. 
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13 


8.50 ±0.21 


0.99 ± 0.28 


0.47 ±0.17 


4th Nearest - highest 50 % 


16 


8.44 ±0.13 


1.16 ±0.32 


0.47 ±0.14 


Number Density - lowest 50 % 


14 


8.37 ±0.15 


1.19 ±0.33 


0.43 ±0.15 


Number Density - highest 50 % 


16 


8.52 ±0.15 


1.02 ±0.22 


0.48 ±0.15 


Halo Mass - lowest 50 % 


13 


8.25 ±0.13 


1.06 ±0.31 


0.42 ±0.16 


Halo Mass - highest 50 % 


16 


8.67 ±0.13 


1.00 ±0.19 


0.36 ±0.12 


Satellites 


16 


8.63 ±0.13 


1.02 ±0.29 


0.46 ±0.15 


Centrals 


18 


8.05 ±0.18 


1.50 ±0.25 


0.30 ±0.14 


M. - Lbui 


Ige.v relation 






All Galaxies 


43 


9.24 ±0.10 


1.11 ±0.17 


0.53 ± 0.07 


4th Nearest - lowest 50 % 


13 


9.17 ±0.14 


1.35 ±0.29 


0.56 ±0.11 


4th Nearest - highest 50 % 


21 


9.34 ± 0.20 


1.09 ±0.25 


0.55 ±0.11 


Number Density - lowest 50 % 


14 


9.30 ±0.17 


0.98 ± 0.24 


0.53 ±0.12 


Number Density - highest 50 % 


20 


9.20 ±0.15 


1.28 ±0.25 


0.59 ±0.11 


Halo Mass - lowest 50 % 


13 


9.17 ±0.18 


1.28 ±0.27 


0.50 ±0.12 


Halo Mass - highest 50 % 


20 


9.31 ±0.16 


0.94 ± 0.24 


0.62 ±0.12 


Satellites 


18 


9.35 ± 0.29 


1.07 ±0.34 


0.63 ±0.15 


Centrals 


25 


9.19 ±0.11 


1.46 ±0.24 


0.49 ± 0.09 



Table 1. The results of linear regression of the black hole scaling relations for environmentally defined subsamples of galaxies. The scahng relations are given 
by Equation[T] and are presented for M.-cr (the full sample and after restricting to only early type galaxies), M,-Mv and M.-M^uigc. t is an estimate of the 
intrins ic scatter in the relation, while a and /3 are the intercept and slope, respectively. The regressions are performed with the liimix.err method of lKellvl 
<2007l) . 
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Figure 3. M. - cr relation of environmentally defined samples of galaxies. 
The sample is divided in two, based on the 4th nearest neighbour distance 
(top), number of galaxies within 5 Mpc and ± 500 kms^^ {second from 
top), the group virial mass (third from top) and whether the galaxy is a 
central or satellite (bottom). The first three indicators are divided by the 
midpoint in the distribution. In these cases, galaxies in dense environments 
are blue and those in isolated environments are red. The lines show the 
linear regression fits of each subsample. 



Sample Ng^,; o /3 e 

Satellites 27 8.34 ± 0.09 4.61 ± 0.58 0.41 ± 0.08 
Centrals 34 8.24 ± 0.07 6.14 ±0.51 0.37 ± 0.06 



Table 2. The results of linear regression of the black hole scaling relations 
for central and satellite galaxies after removing the known massive BCG 
galaxies. The table lists the number of galaxies in each sample as well as 
the intercept (a), slope (P), and intrinsic scatter (e). 



3.1.1 Early types only 

Because of the known morphology-density relation, and the claims 
that the scaling relations are different for different morphology 
classes, it is possible the observed environmental dependence of 
the scaling relations is created by the underlying morphology dif- 
ferences and the change in relative fractions between environments. 
Consequ ently, we restrict only to those galaxies classified as early 
types bv lMcConnell &Mal j2013l) . This reduces the sample to 52 
total galaxies. We redo the linear regressions after redividing the 
sample by the distributions of their environmental indicators. The 
results of these regressions are shown in the second section of Table 
[T]and are plotted in Figure |4] 

Most interestingly, the relations for 4th nearest neighbor, 
group mass and central/satellite give no qualitative change. Indeed, 
both the intercept, a, and the slope, /3, of the relations for the central 
and satellite subsamples are nearly identical whether all galaxies 
are considered or only the early types. Clearly, the environmental 
dependence of the M, - cr relation does not appear to be coming 
simply from the changing mixture of early type and late type with 
environment. However, this environmental dependence may still re- 
flect finer subsamp les of early and late type galaxies. In particular, 
iGrahaml JioOSallb l) have shown that baiTed galaxies have consider- 
able offset from the M. - a relation of non-barred galaxies. 

The only significant change between the full sample and the 
early type sample is the increasing /3 in the case of high number 
density, accompanied by a reduction in [3 with low number den- 
sity. Nonetheless, the early type value of each slope is within the 
uncertainty of the slope determined from the full sample. 



3.1.2 The role of BCGs 

Until now we have treated all central galaxies the same way. 
However, there is some reason to expect that central galaxies of 
massive clusters have fundamentally different formation mecha- 
nisms from normal, isolated central galaxies. To test to what ex- 
tant these BCGs are creating the different scaling relations, we 
will remove them from the sample. The BCGs, as defined by 
iMcConnell & Md l l20I3l) . are the BCGs of Abell 1836 and 3565 
as well as NGC1399, NGC3842, NGC4486 (M87), NGC4889, 
NGC6086 and NGC7768. Again, we redivide the sample by the 
distributions of environmental indicators and perform linear regres- 
sions to each scaling relation. The results are shown in Figure |5] 
and are tabulated in Table |2] The intercepts of the regression have 
no significant change after the removal of the BCGs. In contrast, 
the slope of the centrals becomes shallower by ~ 0.3 after the re- 
moval of the BCGs. However, there is still a significant difference 
in the slope of the relations between satellites and centrals after the 
removal of the known BCGs. 
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Figure 4. M. - c relation for environmentally defined samples of early type 
galaxies. The sample is divided in two, based on the 4th nearest neighbor 
distance (top), number of galaxies within 5 Mpc and ± 500 kms^^ {second 
from top), the group viiial mass {third from top) and whether the galaxy is 
a central or satellite {bottom). The first three indicators are divided by the 
midpoint in the distribution. In these cases, galaxies in dense environments 
are blue and those in isolated environments are red. The lines show the 
linear regression fits of each subsample. 




100 

Velocity dispersion (km/s) 



Figure 5. M, - cr relation for central and satellite galaxies after the removal 
of the 8 known BCGs. The blue and red lines show the results of a hnear 
regression on the centrals and satellites respectively. 



3.2 M.- Mbuig^ relation 

The M, - Albuige relation is potentially a powerful indicator of 
the formation history of bulges. In this section, we examine this 
relation for the 34 galaxies in the sample with good dynamical 
measurements of bulge mass. This sample includes only early type 
galaxies. Using the full sample of 34 galaxies, we find relations 
of a = 8. 44 ± 0.09 and /3 = 1-10 ± 0.17, which are similar to the 
results of iMcConnell & Mai i2013i) . Performing a linear regression 
on several envirormientally divided samples using Equation [T]gives 
the results tabulated in the third panel of Table [T] The fewer num- 
ber of galaxies in this sample prohibits finer bins than dividing in 
half. We plot the results, including the determined regression lines, 
in Figure|6] 

Dividing the sample into isolated and dense environments via 
4th nearest neighbor or 5 Mpc number density leads to no statisti- 
cally significant difference in the scaling relation. In contrast, we 
see that the group mass division leads to an increase in a in dense 
environments. The most massive halos (a = 8.67 ± 0.13) tend to 
host BHs which are higher for a given bulge mass than is the case 
for lower mass halos (a = 8.25 ± 0.13). Despite this offset, the 
slopes of the two relations are indistinguishable. 

The central galaxies (/? = 1.50 ± 0.25) have much steeper 
slope than the satellite galaxies (/? = 1.02 ± 0.29) but also have 
a lower intercept (a = 8.05 ± 0.18 vs 8.63 ± 0.13). Unfortunately, 
in addition to the small sample size, the restriction to a sample of 
early type galaxies results in the centrals and satellites populating 
different sections of the M, - Mbuige relation. While most of the 
central galaxies in this plot have bulge masses greater than 10^^ 
M0 there are many satellite galaxies below that mass. This is not 
unexpected, as it is well known that low mass early type galaxies 
are preferentially satellite galaxies, but it makes it difficult to deter- 
mine what is driving the difference in this sample. Restricting to a 
mass matched sample is not possible with the small numbers. 
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Mbuige (LoQioMsun) 

Figure 6. Mbh - Mbuigc relation of environmentally defined samples of 
galaxies. The sample is divided in two based on the 4th nearest neighbor 
distance (top), number of galaxies within 5 Mpc and ± 500 kms^^ (second 
from top), the group viiial mass (third from top) and whether the galaxy is 
a central or satellite (bottom). The first three indicators are divided by the 
midpoint in the distribution. In these cases, galaxies in dense environments 
are blue and those in isolated environments are red. The lines show the 
linear regression fits of each subsample. 



3.3 M.- Lbuigcv relation 

Another often used BH scaling relation is that between the BH 
mass and the luminosity of the bulge. This is similar to using the 
bulge mass, except it includes the variation of mass to light ratios 
(M/L) of the bulges which can encode information about metallic- 
ity, star formation histor y, dust, etc. Again, we use the sample of 
iMcConnell & Mai j2013l) for the bulge luminosities of 43 galaxies 
in the V band. Our regression of the full sample (a = 9 . 24±0 . 10, 
^ = 1.11 ±0.1 7) is consistent with the lMcConnell & M j (201?) re- 
sults. The results for environmentally split samples are shown in 
Figure |7] and tabulated in the bottom panel of Table [T] 

The dense environments (as measured by 4th nearest neighbor 
or number density), the more massive groups, and central galaxies 
all have steeper slopes and lower intercepts than their correspond- 
ing samples. However, before interpreting these results too deeply, 
we note that the difference between the luminosity distributions of 
the environmentally split samples is even more extreme than the 
bulge mass distributions of the central/satellite relation mentioned 
in Wi.2\ Satellite galaxies are more likely to have truncated star for- 
mation, so they are likely to have lower M/L than central galaxies. 
Thus a central and satellite galaxy may have the same bulge mass, 
but recent star formation in the central would make it brighter. If 
we assume that the power law in Equation [T] is a perfect descrip- 
tion of the data, then w e should not be concerned by the different 
distributions. However, iGrahamI ( 1201 2h has shown that the bulge 
mass and bulge luminosity scaling relations are not well described 
by single power laws. Indeed, it seems clear that the much better 
studied velocity dispersion - luminosity relationship is not a single 
power law, which impli es that the M. - a and M. - My. bulge can 
not both be po wer laws ( lBinnevl[l982l : lMatkovic & Guzmanll 20051 : 
iGrahamllioi ij). It is interesting that we do see a significant differ- 
ence between each of the environmentally split samples; however, 
to tell if this relates to the different luminosity distributions and an 
underlying non-power law relation, or an intrinsic environmental 
difference, requires a more extended sample. Given that high lu- 
minosity satellites are unlikely, the best progress can be made by 
examining a sample of early type, low luminosity centrals (if they 
exist). 



4 DISCUSSION 

In each of the three scaling relations we have examined, central 
galaxies have shown a steeper slope than satellite galaxies. The 
satellite galaxies in the M. - Mtuige relation also showed a higher 
intercept than the central galaxies, while the differences in the in- 
tercepts of the other two relations were not statistically significant. 
The difference in the scaling relations between central and satel- 
lite galaxies was consistentally the most significant of any of the 
environmental divisions. Further, this difference was maintened af- 
ter restricting to early type only samples or removing the known 
BCGs. For these reasons, we will assume that the fundamental dif- 
ference in environmentally divided scaling relations is driven by 
the difference between the central and satellite galaxies in the M, 
- a relation. While the discussion below is framed in those terms, 
the focus on this instead of another relation is largely because ve- 
locity dispersion is a direct probe of the potential well and thus the 
needed escape velocity. 
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Figure 7. Mbh - Ly, bulge relation of environmentally defined samples of 
galaxies. The sample is divided in two based on the 4th nearest neighbor 
distance {top), number of galaxies within 5 Mpc and ± 500 kms^^ {second 
from top), the group viiial mass {third from top) and whether the galaxy is 
a central or satellite (bottom). The first three indicators are divided by the 
midpoint in the distribution. In these cases, galaxies in dense environments 
are blue and those in isolated environments are red. The lines show the 
linear regression fits of each subsample. 



4.1 Secular growth of black holes 

Many authors have investigated the origins of the M. - a rela- 
tion analytically by assuming that the M, mass is due to the sec- 
ular growth of a galaxy without continued cosmological accretion. 
Silk & Rees ( 1998) have shown that an AGN driven outflow which 
conserve s its energy w ill result in an M. - a relation of the form 
g- MKindfiooi) pointed out that because of the short Comp- 
ton cooling time in the inner regions of the galaxy, AGN wind 
shocks would quickly cool and the outflow would be 'momentum- 
driven'. He found that the momentum driven wind leads to a rela- 
tion of the form M.oc a^. After incl uding the adjustment for sub- 
Eddington luminosities suggested bv lMurrav et al.l fcoOSi) . the full 
relation is 



fa K_^4 



(2) 



where k is the opacity and F is the fraction of the Eddington lu- 
minosity at which the black hole emits ( ^^" ). To explain the dif- 
ferential slopes of the centrals and satellites requires an additional 
dependence on a for the centrals but not the satellites. The slope of 
satellites is roughly consistent with a Mbh oc cr*. Thus, the slopes 
imply that ^jA- (x where P ^ 2 for centrals and ^ for satel- 
lites. 

The relevant opacity for Equation |2]is the electron scattering 
opacity. For thermal motions significantly less than the rest mass 
energy, T<< lO'"* K, this opacity is independent of photon fre- 
quency and is simply k = , where ae. is the Thomson cross- 
section. It is unlikely that these terms have a sigificant dependence 
on a or on whether the galaxy is a central or a satellite. 

The static 'protogalaxy' models used analytically implicitly 
assume that the black hole can accrete gas at an unencumbered rate 
which leads to no dependence of F on cr or M, . Nonetheless, the 
Bondi-Hoyle accretion rate is oc M^ while the Eddington rate is oc 
M.. Thus, if the black hole growth is accretion limited, a Bondi- 
Hoyle accreting black hole will have a F variation with its mass. 
However, assuming this F variation does not occur, the last possi- 
bility for a (7 dependence is the gas fraction, fg. 



4.1.1 The dependence of gas fraction on g 

Satellite galaxies infalling into groups and clusters experience ram 
pressure stripping of their hot gaseous halo and, perhaps, of the 
cold gas from the inner regions of the galaxy. These satellite 
galaxies are known to be d eficient in gas and are not forming a 
signifi cant amount of sta rs jHavnes et alii 19841: IXermev & Yound 



1989'; 'Bose lli etal.ll 19971 ; iRasmussen et alj|2012l ; 1jaffe et allllOia 
Fabello et al.l l20ia . Thus, the gas fraction of satellite galaxies is 
likely universally low, with no significant dependence on a. In 
contrast, central galaxies have declining gas fractions with increas- 
ing mass jSaintonge et alj201 ll ; lcatinella et alj2oT2h . This implies 
that the fg oc a~'^ where a depends on exactly how gas fractions 
decline with mass. Unfortunately, this is the wrong dependence to 
explain the observed scaling relations. However, the current f g - cr 
relation is not necessarily the relevant quantity if black hole growth 
has been quiescent for some time. The fg at the time of the last ma- 
jor growth phase could be constant with velocity dispersion, and 
the current f g - cr relation is a consequence of the resultant AGN 

outflows. 

Recently, IZubovas & Kind j2012l ) explored the effect of fg 
changes due to gas depletion from star formation and gas replen- 
ishment in clusters. They suggested that the dominant form of gas 
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loss in both dense and isolated environments is the conversion of 
gas to stars. They argued that the ram pressure stripping of galaxies 
in dense environments results only in the removal of gas on the out- 
skirts of the galaxy, which lowers the star formation rate but does 
not significantly change the central BH mass required to drive an 
outflow. They further argued that cooling flows can replenish cold 
gas in some cluster galaxies (especially BCGs). This would actu- 
ally increase the gas fractions in cluster galaxies and thus require 
higher BH masses at a given velocity dispersion. Thus, they expect 
an offset in the relation for cluster galaxies and isolated galaxies, 
but with the same slopeQ It is worth pointing out that we do see 
this offset when dividing the sample by halo mass, which may ac- 
tually help to explain that relation. However, their requirement that 
satellite galaxies have higher gas fractions than centrals of the same 
mass is hard to reconcile with the observations. It appears difficult 
to explain the different slopes in centrals and satellites by using 
secular growth as a guide. 



4.2 Merger driven growth of black holes 

Mergers are an important part of galaxy growth and a direct con- 
sequence of a ACDM universe. For our purposes, the important 
quantity is the merger rate for centrals and satellites as a function 
of velocity dispersion (or some other mass indicator). Observation- 
ally, the merger rate is always difficult to measure because of its 
transient nature and the uncertain timescale of its disruption fea- 
tures. To add to this uncertainty, low mass satellite galaxies natu- 
rally have more close pairs than similar mass central galaxies and 
thus have substantial 'projection effects'. The merger rate can be 
directly measured in the two environments through cosmological 
dark matter simulations which show that satellite-central mergers 
are more likely f or non-dwarf satellit e galaxies than are satellite- 
satellite mergers jAngulo et al. I l2009h . Assuming that baryons do 
not significantly alter this picture, then satellites live a relatively 
quiet life until they are subsumed by the central galaxy. However, 
given the uncertainty, it is worth exploring whether the nature of 
black hole growth in mergers differs between centrals and satel- 
lites. 

There are two processes important for the black hole scaling 
relations in the context of mergers. Firstly, during a merger, the 
central black holes in each of the galaxies can also merge, or pos- 
sibly be ejected from the dynamics. This will create differential 
satellite/central relations only if the frequency of each occurence 
is different. Potentially the higher gas fractions in central galax- 
ies could damp this effect, leading to more growth in centrals. 
Again, however, this would imply central galaxies have a shal- 
lower slope than satellite galaxies, in contradiction to our results. 
Secondly, if mergers induce bursts of star formation they could 
also drive gas towards the central black hole fo r continued growth 
dSanders & Mirabelll996l : lHopkins et al.l2006al) . If this burst phase 
of star formation driven AGN growth is significantly different from 
the regular growth of black holes, then it could lead to a differential 
evolution of the central and satellite relations. 



1 IZubovas & Kind )2012h actually suggest the existence of 4 offset M, - 
cr relations: isolated spirals, cluster spirals, isolated ellipticals and cluster 
ellipticals which result from the combination of higher fg in clusters and 
larger bulge sizes (requiring longer quasar 'on' phase) in ellipticals. 



4.3 The acclerated growth of black holes in low mass 
satellites 

If we assume that the M, - a relation of centrals and satellites 
does not evolve with redshift, then low mass central galaxies must 
grow their black holes when they become satellite galaxies. The 
linear fits for the full BH sample give a satellite black hole mass of 
2.0x 10^ Mq for a (T = 120 kms~^ galaxy while the corresponding 
black hole mass of central galaxies is slightly more than a third of 
that (7.3 X 10^ Mq). The situation switches at higher velocity dis- 
persion, such that central galaxies of cr = 320 kms~^ have black 
hole masses nearly twice those of satellite galaxies (3.8 x lO^M© 
vs 2.1x10^Mq). 

Given that low mass galaxies have considerable gas fractions, 
it may be that the difference between central and satellite relations 
is driven by these galaxies alone. A plausible scenario for these ob- 
served trends is that low mass (and therefore gas rich) galaxies have 
a phase of increased black hole growth during the transformation 
from a star forming central to a passive satellite galaxy. These pas- 
sive satellites eventually merge with the central galaxy in the halo 
causing the relative increase in black hole mass at high velocity 
dispersion. It is important to remember that most of the merger ac- 
cumulated mass occurs f rom minor mergers w ith objects 1/10 the 
size of the main galaxy jLacev & Coi3ll993l) . These many minor 
mergers with over-massive satellite black holes are likely to lead to 
an over-massive central black hole in massive central galaxies. In 
contrast, high mass galaxies which become satellites do not have 
gas available for an additional growth of the black hole and thus 
have black hole masses similar to centrals. 

We can obtain an estimate for how much of the mass of a 
320 kms~^central galaxy must have come from satellites of 1/lOth 
the mass to explain its m assive black ho l e with out any acclerated 
growth within the galaxy. ' Bemardi et al.l ( |2003|) showed that early 
type galaxies follow a relation such that cx R^'^'', where is 
the Sloan z' band luminosity. Assuming that these galaxies are in 
virial equilibrium and that the z band mass to light ratio is constant, 
then cr cx L°'^^. As such, a galaxy which is 10 times more luminous 
has a ~ 4 times higher velocity dispersion. As mentioned, a 320 
kms~^ central galaxy has a black hole twice the value of a satellite 
black hole at the same velocity dispersion, while a satellite galaxy 
with 1/10 the luminosity, cr = 75 kms~^ , would house a black hole 
6 times larger than that of a central of the same velocity dispersion. 
Thus, if only 20% of the mass of a 320 kms~^ central galaxy has 
been made up by the accretion of satellites of I/IOth the mass, the 
mergers of their black holes will naturally produce a central black 
hole mass twice the size of a black hole in a satellite of the same 
mass. This amount of sa tellite accretion is easily accommodated 
by semi-analytic models dOe Lucia & Blaizo3l2007i ). In effect, we 
require only that gas rich, low mass galaxies to have a phase of 
intense black hole growth to produce the relations we see. 



5 CONCLUSIONS 

We have calculated four environmental indicators for a sample of 
galaxies with direct measurements of black hole mass, velocity dis- 
persion, bulge luminosity and bulge mass. By examining how the 
scaling relations depend on these environmental parameters we ar- 
rive at the following conclusions. 

• The strongest dependence on environment of the M.- cr rela- 
tion is seen when a division by satellite or central designation is 
made. While the vertical offsets of the relation of both samples are 
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similar, central galaxies have a significantly steeper slope ( 6.39 ± 
0.50 vs 4.78 ±0.51). 

• The differential in central/satellite slopes is due not to an un- 
derlying difference in their varying fraction of early type and late 
type galaxies, since the slopes are nearly identical after the removal 
of spirals from both samples. The differential central/satellite slope 
remains even after removing the known Brightest Cluster Galaxies 
from the samples. 

• The steeper central galaxy slope also occurs in the M, - 
Mbuigc and M. - Mv, bulge relations. The smaller sample sizes of 
these relations mean the significance of the trends is less in these 
relations. 

• A simple model of secular black hole growth in which the 
black hole mass is coupled to the current gas fraction predicts that 
central galaxies would have shallower slopes than satelUtes, in con- 
trast to our results. 

• We suggest that gas rich, low mass galaxies experience a pe- 
riod of acclerated black hole growth during the transition from cen- 
tral to satellite. These galaxies, and their large black holes, eventu- 
ally merge with massive central galaxies, leading to more massive 
black holes in massive central galaxies. 

The strong environmental dependence of black hole scaling 
relations provides a new test for models of the co-evolution of black 
holes and their galaxies. We suggest the most progress will be made 
by increasing the number of observations of low velocity disper- 
sion, bulge mass and bulge luminosity isolated central galaxies. 
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